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Abstract. Consider a homogenous fluid membrane, or vesicle, described by the 
Hclfrich-Canham energy, quadratic in the mean curvature. When the membrane is 
axially symmetric, this energy can be viewed as an 'action' describing the motion of a 
particle; the contours of equilibrium geometries are identified with particle trajectories. 
A novel Hamiltonian formulation of the problem is presented which exhibits the 
following two features: (i) the second derivatives appearing in the action through the 
mean curvature are accommodated in a natural phase space; (ii) the intrinsic freedom 
associated with the choice of evolution parameter along the contour is preserved. As 
a result, the phase space involves momenta conjugate not only to the particle position 
but also to its velocity, and there are constraints on the phase space variables. This 
formulation provides the groundwork for a field theoretical generalization to arbitrary 
configurations, with the particle replaced by a loop in space. 



PACS numbers: 87.16.Dg, 68.03.Cd, 02.40.Hw 

Phospholipid molecules self- assemble in water to form vesicles, or membranes P]I2|- 
The vesicles are very thin compared to their size so that it is sensible to describe them 
as surfaces of negligible thickness. Moreover, the vesicle behaves as a two-dimensional 
fluid: there is no resistence to shear so that the molecules move freely in the plane of 
the membrane. In a theoretical model of lipid vesicles, this physical property is realized 
as invariance under reparametrizations of the surface. The energy thus depends only 
on the geometry of this surface. The leading term in the energy is proportional to the 
integrated square of the mean curvature, penalizing bending [HIIHE]- The molecular 
details are essentially irrelevant. 

The 'shape equation', describing equilibrium geometries [HllZj is a fourth order non- 
linear PDE. With axial symmetry, the PDE reduces to a non-linear ODE which, in turn, 
possesses a first integral [HI E] • Indeed, it is possible to interpret the energy as an action 
describing the motion of a particle. Geometrical contours can be identified with particle 
trajectories. Building on Deuling and Helfrich's pioneering work in the early seventies 
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[TU] . axisymmetric solutions of the shape equation describing an isolated vesicle were 
pretty well understood by the mid-nineties (see the reviews [TU 112^ IT3]). 

Without the symmetry we are less well off. However, with increasing computer 
power, impressive results can be achieved; Monte-Carlo and dynamical triangulation 
(see e.g. E3) will minimize the energy for us. In the latter case, the program 
Surface Evolver was designed with exactly this sort of problem in mind [16]. 

At this level, the shape equation is consigned to the status of a curiosity. However, 
there is a lot of information encoded in the shape equation which can be accessed 
without having to solve it explicitly. For example, it is not widely known that the 
shape equation can be cast as a conservation law for the stresses prevailing within the 
membrane [Hj. These stresses are completely geometrical. They transmit forces. It 
would be difficult to understand the nature of these forces without taking the shape 
equation apart. In this respect, a computation scheme to solve the shape equation, 
would be a useful complement to energy minimization. The mechanical analogue of 
the axially symmetric shape equation is most naturally formulated as a Hamiltonian 
initial value problem. There is no obstacle, in principle, to setting up a field theoretical 
generalization: instead of a point particle take a closed loop; motion of the loop will 
generate a surface. Unfortunately, the existing Hamiltonian approaches to solving the 
axially symmetric shape equation, that use arc-length along the contour as a parameter, 
are tailored very specifically to the symmetry, so they are not very helpful. 

In this paper, we present a novel Hamiltonian formulation of the axially symmetric 
shape equation which takes no shortcuts home. It will, however, admit a field theoretical 
generalization with the particle replaced by a loop in space JJTj. This formulation will 
involve two key features: 

(i) When axial symmetry is relaxed there is no single priveledged parameter analogous 
to arc-length along the contour. The formalism should therefore respect the intrinsic 
freedom associated with the choice of evolution parameter. 

(ii) A point that tends to go unnoticed in the axially symmetric is that the action 
involves not only first derivatives (velocities) but also second derivatives (accelerations), 
a feature that is somewhat challenging from a Newtonian point of view. With axial 
symmetry, the problem is simply sidestepped by introducing the turning angle along 
the contour ( a velocity) as an intermediate variable; with respect to this variable, the 
action involves no derivative higher than first. What amounts to the same thing, only 
without the sleight of hand, is to introduce the natural phase space that is appropriate 
for the Hamiltonian formulation of a theory based on an action involving accelerations: 
introduce momenta not only canonically conjugate to the particle position, but also to 
its velocity. 

Even if axial symmetry were to be our final goal, there are benefits to this apparently 
unnecessarily complicated formalism: both the momenta and the constraints possess 
physical meaning and Hamilton's equations will evolve physical initial data in a 
remarkably straightforward way. As we will show the difficulty is in the setup; not 
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We model a lipid vesicle as a two-dimensional surface E. The surface is described 
locally by the embedding x = X(w a ), where x are local coordinates in space, u a = 
(u l ,u 2 ) local coordinates on the surface, and the position functions X(w a ) are three 
functions of two variables. We denote by e a = <9 a X = dX/du a the two tangent vectors 
to the surface. The metric induced on £ is given by their inner product, g ab = e a ■ e^. 
The unit normal n to S is defined implicitly by e a • n = 0, n 2 = 1. The extrinsic 
curvature tensor is K a j, = — n ■ <9 a <9feX, and the mean curvature is K = g^K^, where g ab 
is the inverse of the induced metric g^. In terms of the principal curvatures, {c\, c 2 }, we 
have K = c± + C2- The intrinsic scalar curvature can be given in terms of the extrinsic 
curvature via the Gauss-Codazzi equation as 7Z = K 2 — K ab K ab ; it is twice the Gaussian 
curvature G, i.e. 7Z = 2G = 2c\C2- 

We consider the Helfrich-Canham geometric model, or bilayer coupling model, for 
a fluid lipid vesicle, with energy 

F[X] = ^JdAK 2 + (3jdAK + aA-PV. (1) 

where the constant Hi is the bending rigidity, dA = ^fgd 2 u denotes the infinitesimal 
area element on the surface, and g is the determinant of the induce metric g a b- The 
constants /3, cr, P are Lagrange multipliers enforcing the constraints of constant total 
mean curvature (or constant area difference between the layers), constant area and 
constant enclosed volume V, respectively ^H]- A refinement, known as the ADE model, 
imposes a non-local constraint involving the square of the area difference JTHJ 1201 EI] 
Our considerations can be extended to this and other geometrical models for membranes. 
Note that the volume can be written as a surface integral: 

V = - f dA n • X . (2) 



3 

We have not included a term corresponding to the Gaussian bending, Pb[X] = 
k g J dATZ, since it is a topological invariant by the Gauss-Bonnet theorem, and, as such, 
it does not contribute to the determination of equilibrium configurations. The energy ([!} 
is invariant under rigid motions, translations and rotations, of the surface in the ambient 
space. It also possesses a local symmetry: invariance under reparametrizations. 

The vanishing of the first variation of the energy (0), with respect to variations of 
the position functions X(w a ) — > X(n a ) + 5X(u a ), gives the shape equation jBJ 1221 

A " ' ~ 1 + (371 + oK - P = , (3) 



Hi 



-V 2 K - y(K 2 -271) 



where V 2 denotes the surface Laplacian. This fourth order non-linear PDE determines 
the equilibrium configurations of lipid vesicles. There is only one equilibrium condition, 
whereas naively one would have expected three. Reparametrization invariance informs 
us that two linear combinations of these three equations, corresponding to tangential 
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deformations, must vanish identically [2S1- The only physical deformations are those 
normal to the surface. 

Let us now specialize to axially symmetric configurations. The embedding of an 
axially symmetric configuration can be written as 

/ R(t) cos0, \ 

x = X(w a ) = X(t,0) = R(t) sin0, , (4) 

V Z(t) ) 

where t is an arbitrary parameter along the contour of the surface at fixed 0. Any space 
vector V can be written in adapted components as 

/ V R (t) cos0, \ 

V(t,0)= V R (t) sin0, , (5) 
V V z (t) j 

so that on the plane = it reduces to a two-dimensional t-dependent vector with 
components {Vn(t), Vz(t)}. The basis adapted to the surface is given by the two tangent 
vectors e t = d~K/dt = X and = <9X/<90, together with the unit normal vector 



/ Z(t) cos ( 
n(t,<f>) = — Z(t) sin< 

V -m 

where we introduce the function 

N = 



\ 



R 2 + Z 2 . 

Note that arclength I along the contour is defined infinitesimally by dl 
induced metric and the extrinsic curvature tensor assume the form 



(6) 



(7) 

N dt. The 



9ab 



N 2 






R 2 



K, 



ah 



1 ( RZ -ZR 
N [ RZ 



For the mean curvature and the scalar curvature it follows that 
R(RZ - ZR) + N 2 Z 2Z(RZ - ZR) 



K 



g ab K ab 



n 



(8) 



(9) 



RN 3 ' " RN 4 

The Helfrich-Canham energy (0) specialized to a axially symmetric configurations, 
in terms of an arbitrary parameter t, is 

(10) 



F[X] = 2tt / dt L(R, Z, R, Z, R, Z) , 
where the Lagrangian function is 



L(R, Z, R, Z , R, Z) 



k [R(RZ - ZR) + ZN 
2 RN~ 5 



212 



R 



+ (3j^(RZ-ZR) 



(3Z 



+ aRN - -R{RZ - ZR) . 



;n) 



We now treat this energy as an action determining the motion of a fictitious particle in 
the two dimensional configuration space {R, Z}. The arbitrary parameter t will play the 
role of time. The Lagrangian function involves the squared acceleration of this fictitious 
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particle; it enters quadratically due to bending, and linearly due to the constraint of 
constant total mean curvature. The only serious non-linearity is the dependence of 
the Lagrangian on the velocity {R,Z}. The factors of R,N ensure that the action 
is invariant under reparametrizations of t, all that remains of the reparametrization 
invariance of the energy (0) once we specialize to axially symmetry. The only dependence 
on Z is through the volume, if P 7^ 0. 

The axially symmetric version of the shape equation can be obtained either by 
direct specialization of the general shape equation 0, or as the vanishing of the Euler- 
Lagrange derivative of the energy ()1U|) (see the Appendix) 

F _^_(dL\ d_(dL\ dL_ 

£l ~^\^)~Jt\dQ^) + W ) ( J 

where Q % = {R, Z}. Each gives a copy of the shape equation. 

Two features of the dynamical system defined by the action complicate its 
Hamiltonian formulation: the energy depends on second derivatives of the configuration 
variables, the position functions; and because we have chosen to work in an arbitrary 
parametrization, there is a local symmetry - reparametrization invariance. The 
implications of the latter will become apparent only after we have obtained the canonical 
Hamiltonian. First we must construct the phase space. 

The most direct approach to handling the presence of the second derivatives {R, Z} 
is to extend the phase space: we treat not only the position variables {R, Z} but also 
their velocities {R, Z} as configuration variables, and introduce conjugate momenta 
for both sets of variables. (A brief summary of the Hamiltonian formulation of higher 
derivative systems is provided in the Appendix.) The momenta {Pr, Pz} conjugate to 
the velocities {R, Z} are, respectively, 

F) T 7 R 7 

PR= m = - K w [R(RZ ~ ZR) + ZN2] ~ ' (13) 

Pz = ^ = n^[R(RZ - ZR) + ZN 2 ] + . (14) 

We note that the vector {Pr, P z } is directed along the normal to the contour, and that 
its bending part is proportional to the mean curvature [T7]. Note also that the vector 
{Pr, Pz} depends at most on second derivatives of {R, Z}. 

The momenta {pr,Pz} conjugate to {R, Z} are, respectively, 

dL _ d (dL\ 
PR ~dR~~dt \dRJ 

5kR[R(RZ - ZR) + ZN 2 ] 2 k(RZ + 2RZ) [R(RZ - ZR) + ZN 2 ] 



2RN 7 RN 5 
2(3RR(RZ-ZR) (3RZ aRR P ov . 
m + 1^ + — + 3 RZ - Pr > (15) 

PZ ~dZ~di [dZ 
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5kZ[R(RZ - ZR) + ZN 2 } 2 + k(-RR + R 2 + 3Z 2 )[R(RZ - ZR) + ZN 2 



2RN 7 RN 5 
23RZ(RZ - ZR) 3RR o RZ P 2 • 

m "A^ + ^r + /3 ~3^ ~ Pz - (16) 

Despite the unpromising appearance of these expressions, in J7| we will see that the 
vector {pr,Pz} is the projection of the stress tensor associated with the membrane along 
the unit tangent to the contour. There is a direct physical significance attached. Note 
that this vector has a dependence, through the derivatives of {Pr,Pz}, on the third 
derivatives of {R, Z}. 

We have now identified the appropriate phase space for the system defined by the 
energy (fTU|) : the position of a particle in two dimensions {R, Z}, and its conjugate 
momenta {pr,Pz}, given by (fT^j) . (fTtjj) . together with the velocity {R, Z} and its 
conjugate momenta {Pr,Pz}, given by |T3|) . |T4j) . Intuitively, the position is conjugate 
to its third derivative; the velocity is conjugate to the second derivative of the position. 
As we will see below, however, not all of this phase space is accessible: reparametrization 
invariance will imply constraints. 

Our next step is to construct the Hamiltonian on the phase space. Because the 
Lagrangian depends on second derivatives of {R, Z}, the definition of the canonical 
Hamiltonian H involves the Legendre transformation with respect to the accelerations 
{R, Z} as well as the velocities, {R, Z}, (see the Appendix) as 

H {P R , p R , R, R- P z , Pz , Z, Z) = P R R + P Z Z + p R R + p z Z - L(R, Z, R, Z, R, Z) . (17) 

The definition of the momenta {Pr, Pz} is used to express the higher derivatives 
configuration variables {R, Z} in terms of the phase space variables Pr, Pz, R, Z, Z, R. 
Unlike a lower order Hamiltonian system, the terms PrR, pzZ are left alone; they are 
already in canonical form. 

To facilitate the elimination of {R, Z} in (fT7j) . we square (JTHJ) and ((Hj) , defining 
the momenta {Pr, Pz}, and add to give 



7V 8 

R 2 {RZ - ZR) 2 = 



K 



Z 2 RZ\ f ZR RR\ 

Pr + n— + 3 — — + [ P z - k — - - 3 — - 

, N 3 N 2 \ N 3 N 2 I 



(18) 



It follows that the canonical Hamiltonian is expressed in terms of the phase space 
variables as 



Ho =p R R + p z Z + 



2kR 



2 




hZ 3Z — aRN + —R(RZ — ZR) . (19) 



2RN ' 3 
This Hamiltonian is quadratic in {Pr, Pz} and linear in {pr,Pz}- 

We have dealt with the first difficulty, the higher order nature of the system 
defined by (jl()|) . now we face the second one, the presence of a local symmetry. In 
this higher derivative model, the presence of reparametrization invariance implies that 
the Hessian of the Lagrangian with respect to the second derivatives is degenerate, and 
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it is impossible to invert for the accelerations in terms of their conjugate momenta. The 
Hessian is 

_&L_ kR(Z* -RZ\ 
13 <)()'()' N 5 { -RZ If ) ' 1 ' 

with Q l = {R,Z}, and we see that its determinant vanishes. This means that at any 
value of the parameter t the phase space variables are not all independent, they are 
connected by constraints. The first (or primary) constraint is easily identified from the 
definition of the higher momenta (|13|) . (J 14)) as 

C = P R R + P Z Z = . (21) 

This is simply the statement that the vector {Pr,Pz}, is directed along the normal to 
the contour. As we will, see this is equivalent to the fact that the tangential component 
of the acceleration is gauge: the parametrization we choose will fix this component. 

The Hamiltonian that generates the motion is given by adding this constraint to 
the canonical Hamiltonian, 

H = H + XC- (22) 



the Lagrange multiplier A is an arbitrary function of t that enforces the constraint (j21j) . 

The Poisson bracket appropriate for this higher derivative model is, for any two 
phase space functions /, g (see the Appendix) 

r , i df dg df dg df dg df dg . , . . 

dRoP R ORdpR dZoP z dZ dp z 
the time derivative of a phase space function / is given by the Poisson bracket with the 
Hamiltonian 



f = {f,H} = {f,H } + X{f 1 C}. (24) 

We have identified a constraint C on the phase space variables. This is not the 
whole story, however. Even if C = initially, we are not guaranteed that it continues to 
hold. Consistency requires that C = be preserved by the evolution: a short calculation 
gives 

C = {C,H } = -H . (25) 

Thus we need to impose the secondary constraint 

#o = 0; (26) 

the canonical Hamiltonian itself must vanish, the hallmark of reparametrization 
invariance. Here, it shows up as a secondary constraint. Clearly Hq = 0. There 
are no other (tertiary) constraints. As a constraint, H = specifies the tangential part 
of the vector {pr,Pz} in terms of the remaining dynamical variables. 

The Hamiltonian function ()2*2*|) generating the dynamics is a linear combination of 
two constraints. Hamilton's equations will reproduce the equilibrium condition given 
by the vanishing of the Euler-Lagrange derivative p2jl . 
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The first pair of equations is 
dR dH 



R, (27) 



dt dpR 
dZ dH 

n = wr z - (28) 

since {pr,Pz} appear in the Hamiltonian only in the combination PrR + pzZ. These 
equations tell us how the vector {R, Z} evolves; in this formalism they are model 
independent. 

The second pair of equations is 

dR ~ dH iV 3 / Z 2 n RZ\ • 

^- R= ap-r^{ P " + K ^ + ^) +XR ' (29) 

dZ - dH N 3 I n ZR „RR\ 

They tell us how {R, Z} evolves. They involve the Lagrange multiplier A explicitly. 

Just as (j2*Tj) and (J28j) encode the definition of the canonical variables {R, Z} as the 
time derivatives of {R, Z}, one would expect these equations to encode the definition of 
the momenta {Pr, Pz} in terms of {R, Z}, {R, Z} and {R, Z}. 

Let us first express the Lagrange multiplier A in terms of the acceleration. We 
multiply (j2Uj) by R and (J5U|) by Z and we add. Using the constraint (j2*Tj) . we identify 

_ RR + ZZ N 
X ~ N 2 ~2N*- (3) 
It vanishes in a parametrization by arc-length. Geometrically, it is the affine connection 
for the planar curve described by {R(t), Z(t)}; the component of the acceleration tangent 
to the contour is pure gauge-it can be chosen arbitrarily. In particular, it can be chosen 
to vanish. If the expression (|3*Tj) for A is fed back into (|2*H|) and (|3()jl. we find that they 
reproduce the form (fT3~j) and (fTIj) for P R and Pz, respectively. We do, however, have to 
use the primary constraint. 

The third pair of equations is 

dP R dH xn 3fflV /n2 n2 x P Z Z f3 2 RR 



PR 



\P R - — (P 2 R + P 2 Z 
H 9kP \ R z 



dt OR 2kR \ k *J R 2kN 

+ -^[Pz(2R 2 + Z 2 ) - P R RZ] + ^ + ^RZ , (32) 
dPz dH xn 3ZN f 2 <a P Z R 2P R Z (3 2 RZ 



dt dZ 2kR v " *J R R 2kN 

- -^l P n(R 2 + 2 ^ 2 ) " P zRZ] + ^ ■ (33) 

They tell us how the vector {Pr, Pz} evolves. One would expect these equations to 
encode the definition of the momenta pr and pz given by (JT5j) and (JTBj). To show this 
is not entirely straigtforward. It is necessary to use the information gathered in the 
previous Hamilton equations, namely the form of {Pr, Pz} and of A. 
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Finally, the fourth pair of equations is 

dpn _dH N 3 P R Z 2 P z RZ (3*N 

dt dR 2kR 2{ r z) R 2 R 2 2k 

+ crN -^(2RZ - ZR) , (34) 

dp z _ dH_P ■ 

-df--dZ-3 RR - (35) 
They tell us how {pr,Pz} evolve. With these equations, we reproduce the vanishing of 
the Euler-Lagrange derivative (fT2j) . One sees that the latter of the two equations has 
the obvious first integral 

p 

J = p z - -R 2 = const. (36) 
6 

The first integral of the axially symmetric shape equation [HJ H| emerges naturally within 
this framework. 

The recipe to construct an axially symmetric equilibrium configuration is as follow: 
Choose initial data: At t — 0, take a point on the plane, specified by its position vector 
{R, Z}, choose a velocity vector {R, Z} (this encodes the initial direction of the contour); 
next choose a vector {Pr, Pz}, orthogonal to the velocity (so as to satisfy the primary 
constraint ipTjl ): finally choose the momentum {pr,Pz} with a tangential component 
consistent with the secondary constraint Hq = where Hq is given by Eq. (fT9"|) . These 
are our physical degrees of freedom. 

This initial data set is evolved using Hamilton's equations. An equilibrium surface 
contour {R(t), Z(t)} will be generated. The contour itself will not depend on the choice 
of the lagrange multiplier A (or equivalently the choice of the parameter t). 



In this paper we have examined the construction of axially symmetric equilibrium 
configurations of a fluid membrane described by the Helfrich-Canham energy from a 
Hamiltonian point of view. If axial-symmetry were our final aim this would be a very 
heavy-handed approach to the problem. The value of all of this formalism will become 
apparent when we consider the generalization to non-axially symmetric configurations 
[TT] . This will involve stepping up from the Hamiltonian dynamics of a particle to the 
corresponding dynamics of a field describing a closed curve in space. The membrane 
surface will be generated by the evolution of this curve. 
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APPENDIX 

For the benefit of the reader unfamiliar with the Hamiltonian formulation of higher 
derivative systems, we consider in this appendix the Hamiltonian description of a toy 
model: a particle moving in one-dimension described by a Lagrangian of the form 
L = L(q, q, q). The Euler-Lagrange derivative for this Lagrangian is 

£= *_(dL\_ d_ (d_L\ d_L 

dt 2 \dq) dt\dq) dq' 1 ' 

The phase space is given by the two conjugate pairs {q, P} and {q,p}- The momenta 
conjugate to q and q are, respectively, 

dL d (dL\ , . 

(39) 

The canonical Hamiltonian is constructed as the Legendre transformation with respect 
to both the acceleration q and the velocity q as 

H(q,P;q,p)=Pq + pq-L, (40) 

where one uses the definition of the higher momentum P to express the highest derivative 
q in terms of the phase space variables P, q, and q. The term pq in the canonical 
Hamiltonian is left alone, since it is already in canonical form. 

The Poisson bracket appropriate for this higher derivative model, for two arbitrary 
phase space functions /, g, is 

ffX dfdgdfdg 

{L9}= d-qdP + d-qdp-- {f ~ 9) > (41) 
and the time derivative of a phase space function is given by this Poisson bracket with 
the Hamiltonian 

! = {/>#}■ (42) 
In particular, it follows that the Hamilton equations are 
dq OH . 

Tt=^p- = q > (43) 

d A = d JL (44 ) 

dt OP ' 1 ] 

* = ~w (45) 

? = -^- (46) 
dt dq y ' 

The first equation identifies the time derivative of q with the canonical variable q; the 
second equation identifies the form of the momenta P conjugate to q; the third equation 
identifies the momenta p conjugate to q modulo the definition of P. Using the first three 
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equations, the fourth equation then reproduces the vanishing of the Euler-Lagrange 
derivative (pT7J). 

One important special case is given by a Lagrangian linear in the acceleration, 
L = g(q, q)q. In this case, the higher momentum P = g(q, q) is independent of q so that 
the acceleration cannot be expressed in terms of the canonical variables. However, it is 
always possible to add a total derivative to the Lagrangian, and obtain a Lagrangian 
that depends at most on q (see e.g. 
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